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The classical cloud top entrainment instability condition is 
commonly formulated as 

g A 8 e + h Aq^ < 0 (1) 

where g,h are slowly varying parameters taken to be constant, 

&0 e , Aq^ the jumps in equivalent potential temperature and total 
water mixing ratio at the cloud top. 

Plotted in 8 e , q t coordinates the entrainment instability 
criterion (El) for a collection of different atmospheric 
conditions appears as a highly elliptical point cluster as a 
result of the strong correlation between 0 e and q t . More 
independent coordinates are 6% and q^-, where 8% is the liquid 
water potential temperature or saturation point potential 
temperature. In this paper we examine El in 8^ , q-^ coordinates. 
Rather than focussing on jump conditions we will emphasize the 
gradients Sd^/Sq^ . El occurs when 

t < (& ^ H /& 9t) crit 

where the critical value represents the gradient along the moist 
virtual adiabat. (2) is entirely consistent with (1), however the 
manner in which El is shown by (2) reveals that additional 
parameters need to be considered in the entrainment of negatively 
buoyant parcels. For simplicity we neglect radiation. 

Figure 1 shows a 8 it q t plot typical of 6/29/87 (C130-data) . M 
corresponds to the in-cloud parcel, T to the parcel above the 
cloud. Cloud top appears at 870 hPa. Mixtures of T and M fall 
along the straight line TM. This diagram can be interpreted in 
pressure (p) and saturation point pressure (p*) coordinates. The 
p=870 hPa line shows the mixing level and whether M, or T or both 
are cloudy or clear. Other mixing levels can be arbitrarily 
defined by simply shifting the p-line up or down. If a point, say 
M falls below that p-line it means that at that pressure level a 
parcel with the coordinates M is cloudy. If it falls above that 
line, it is clear. The p*— lines indicate how much lifting or 
lowering is needed to make M or T just saturated. For example if 
the cloudy boundary layer would be well mixed, it would be 
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represented by M and the cloud thickness would be (937-870)= 67 
h p a. In Figure 1 the moist virtual adiabat falls to the right of 
the mixing line MT. This means that all cloudy mixtures (below 

the level p=870 hPa) are denser than the unmixed cloudy parcel M. 

Additional information about the mixing process is available if 
we draw dry virtual adiabats on this plot as in Figure 2. The 
moist virtual adiabat through M kinks at the p=870 hPa line at 
point c and then follows the dry adiabat of 0 VU =291.72 K. T lies 
on the dry virtual adiabat marked as 299.27 K. Therefore there is 
a difference of A0 V = 7.55 K between the dry air and the cloud 

top so that the cloud is statically very stable. In order to make 

the cloud statically unstable it would have to be raised several 
hundred hPa / s until point c would be above the dry virtual 
adiabat through T. Entrainment at pressure level p=870 hPa takes 
parcels with p*=566 hPa at T and mixes them with parcels at p* = 
937 hPa. During the mixing process the saturation point 
coordinates of the mixture slide from T along the mixing line MT 
towards M. Entrainment will increase p* of the mixture until it 
gains the same buoyancy as the cloudy parcel M at point a, where 
the dry virtual adiabat 0 VU =291.72 K intersects the mixing line. 
At a the parcel is still clear. Subsequent mixing from a to b 
decreases the mixture buoyancy below that of the cloudy parcel. 

At b where P*mixtuj-e = Pm;j.xing level the mixture becomes cloudy. 
Below on the mixing line, lines of equal buoyancy are 
represented by lines parallel to the moist virtual adiabat. In 
continuing the mixing process from b to M the buoyancy will 
increase again until it becomes the same as that of the cloudy 
parcel M. 1 

Four points can be immediately recognized from this plot: First 
the densesst mixture is just cloudy (point b) . Although this is' 
well known from the literature, the actual buoyancy difference 
can be immediately read from this graph by drawing the dry 
virtual, adiabat through point b (0 VU =291.45 K) and computing the 
difference with 0 VU through a and c (0.27 K) ; Second, the 
fraction of dense buoyant mixtures to all possible mixtures is 
given by (a-T)/(T-M); Third, the buoyancy difference of cloudy 
mixtures with M can be increased if the unmixed dry parcel would 
be cooler or dryer than represented by point T; Fourth, the 
fraction of buoyant mixtures can be increased if the mixing level 
is raised (pressure lowered) . Lidar cloud top data shows that 

® flight leg cloud top variability of several hundred 

meters (several tens of hPa's) is not uncommon. The two dashed 
lines parallel to the p=870 hPa line indicate how variations in 
the mixing level changes the fraction of potentially positivelv 
buoyant parcels. * 


Neglecting radiation we can derive the virtual potential 
temperature flux from this diagram as follows: Let F,™ be the 
convective flux in (0 e , 6 ^ , q^) and w e be the entrainment 
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velocity, then 

f MT =(T-M)w e =(T-b)w e +(b-M)w e (3) 

where T = ( 6 e , 6 q^- ) ij. 

M = ( 9 e ' 9 A ' St ) M 

then the virtual potential temperature flux is defined as 

F 0v = S u ( 9 AT ~ 9 Ah ) u e + h u ( StT “ Stb ) w e + 

g c ( 9 Ah ~ 9 AM ) w e + h c ( Stb “ StM ) w e ( 4 > 

where g u , hy, g c , h c are the appropriate factors for clear and 

cloudy conditions. From (4) it is immediately clear that although 
a portion of the buoyancy flux ( in cloud ) is indeed positive by 
virtue of the slope of the mixing line MT, this flux is very 
small in comparison to the flux necessary to create the cloudy 
mixture in the first place. We believe that this is a point that 
is commonly overlooked in El-studies. The energy necessary to 
create the cloudy mixture is very much greater (proportional to 
the difference (299.27-291.45 K) than the energy created by 
mixing inside the cloud (proportional to the difference 
291.72-291.45 K) . It is therefore not surprising that recent 
studies have shown that in many cases where condition (1) or (2) 
was satisfied, the clouds appeared to be stable. 

So far we only considered mixing at pressure level p. However a 
dry parcel from above the cloud is drawn into the cloud and 
acquires a downward speed representative of the in— cloud 
circulation. Let this speed be indicated by wip. wip is responsible 
for lowering the mixing level of the parcel, increasing the 
pressure level. On the other hand w e , the entrainment velocity is 
responsible for increasing the saturation point pressure of the 
mixture. The process of entrainment and vertical movement of the 
parcel is schematically represented in Figure 3 . Assume that the 
highest cloud tops are at point b, and that at that level the dry 
unmixed parcels have their saturation level at T (566 hPa) . As 
the parcel is drawn into the layer the mixing pressure level is 
increased from point b to say point a. If point a corresponds to 
the lowest cloud top we know that at point a all mixtures are 
cloudy. This means that in moving an unmixed parcel at point b to 
point a the saturation level of the mixture has increased from 
its unmixed value at T to the actual pressure level at a. Simple 
geometry on Figure 3 reveals the following constraint: 

f = w T / w e < ( p a - pj-> ) / ( p a - Pt ) ( 5 ) 


For typical cloud top variations of 15 hPa we find f < 0.05 
It means that the entrainment velocity needs to be very much 
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higher than the translation velocity in order to get dense 
sinking cloud parcels. This statement means that cloud top 
entrainment instability is unlikely to break up a cloud in case 
the circulation speed near the cloud top is large; in such cases 
the strength of the circulation merely draws unmixed parcels into 
the layer that do not have the chance to become cloudy in their 
downward transport. Below the mean cloud base those parcels will 
always be less dense than the mixed layer environment. 

In conclusion we have shown that a saturation point diagram can 
be used to investigate the details of mixing in cases where the 
cloud top entrainment instability criterion is satisfied. We 
find that for typical situations found during FIRE, where the El 
condition is satisfied clouds are likely to be stable because 
energy required to create a cloudy mixture is much greater than 
the energy which is released once the parcel has become cloudy. 
The mixing level is crucially important in determining the 
fraction of cool cloudy mixtures. The vertical levels at which 
cloud tops can be found (derived from cloud top lidar data) puts 
a constraint on the entrainment velocity and the transport 
velocity of the mixture in such as way that the entrainment 
velocity needs to be an order of magnitude higher than the 
transport velocity in order to get any cloudy sinking mixtures. 

If the transport velocity is too large, clouds will break up, 
however mixtures will remain less dense than their environment. 
Radiation was neglected in this analysis, but is likeky to 
enhance the instability, as it cools the parcels once they become 
cloudy. A complete analysis then involves another velocity scale 
representing the speed at which a mixture moves along lines of 
equal 9 ^ . The mixing process will deviate from the simple mixing 
line structure as shown in this paper and is the subject of 
further research. 
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